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OPTIMUM DESIGN CONSIDERATIONS OF A GUST 


ALLEVIATOR FOR AIRCRAFT 

Waldo I. Oehman 
Langley Research Center 

SUMMARY 

Optimum design considerations of a gust alleviation system for aircraft flying in 
turbulent air are presented in this paper. A vane sensor (with noise) was used to mea- 
sure vertical gusts, and elevators and flaps were used to reduce the root-mean-square 
(rms) value of the normal accelerations associated with the aircraft response to gusts. 
Since turbulence has stochastic properties, stochastic control theory was used in the 
analysis. A performance -index function involving normal acceleration and control 
deflections with a weighting matrix was used. 

Application of the analysis was illustrated by a short take-off and landing (STOL) 
airplane in flight through turbulent air. Effects of varying the noise characteristics of 
the vane sensor and of the weighting matrix in the performance -index function were 
determined. 

Stochastic control theory was applied to calculate the rms response of the airplane 
to turbulence. The calculations showed that a weighting number appearing in the perfor- 
mance index should be less than 10. A value of 3 was selected for subsequent calculations. 
Normal acceleration could be reduced by 92 percent when the intensity of the measure- 
ment noise had a very small value. However, the filter gains were large, and large con- 
trol deflection angles were required. Good alleviation was calculated when the intensity 
of the measurement noise was about 3.6 percent of the vane deflection angle. Normal 
acceleration was reduced by 63 percent for moderate values of the gains. 

INTRODUCTION 

One of the purposes of gust alleviation systems for aircraft is the reduction of nor- 
mal accelerations caused by gusts. Analyses of such systems have made use of various 
mathematical techniques to represent the gust spectra and to evaluate system response. 
For example, reference 1 describes atmospheric turbulence as sinusoids having unit 
amplitude and varying frequencies. The system was analyzed by using the amplitude 
and phase angle of the airplane response to the turbulence as a function of gust frequency. 
In references 2 and 3, atmospheric turbulence was assumed to have a von Karman power 



spectral density function. Stochastic control theory was used to evaluate the root -mean- 
square (rms) responses of the airplane. 

Since a turbulent atmosphere is stochastic , the application of stochastic control 
theory to gust alleviation problems appears to be natural. Reference 4 gives the first 
results obtained in the application of stochastic control theory to flight control problems 
and emphasizes the problem of stability derivative identification during flight in turbu- 
lence. Results concerning gust alleviation were primarily qualitative. 

In this study, stochastic control theory (see ref. 5) is applied directly to the opti- 
mum design of a gust alleviation system for aircraft. Atmospheric turbulence, in this 
study, is assumed to be a random process characterized by a Dryden power spectral 
density function. An angle -of -attack vane was mounted ahead of the wing. The deflec- 
tion angle of the vane was measured and the resulting signal, which had superimposed 
noise, was used to actuate a wing flap and the elevator. The vane signal and the control 
deflections were the inputs to a Kalman-Bucy filter. This filter gives the best estimate 
of the state of the system. Linear combinations of the estimated state variables were 
used to actuate the controls. The performance of the system was evaluated by the per- 
cent reduction of the normal acceleration of the airplane. Quantitative results are pre- 
sented to show the influence of two unknown parameters on the system performance. One 
parameter is a weighting number appearing in the performance index; the second param- 
eter is the intensity of the measurement noise. 

SYMBOLS 


A system matrix, 4x4 

a 2 ;(t) normal acceleration, m/sec^ 

a.Q,a.^,y-^,y 2 parameters used in appendix A 
B input matrix, 4x2 

C measurement matrix, 1x4 

Cjn pitching-moment coefficient, — ;; 

^oo®W^ 

„ „ t ff • 4 . Force in Z -direction 

C 2 Z -force coefficient, 

^oo®W 

c wing mean aerodynamic chord, m 


2 



II nil 


II 


I 


D controlled variable matrix, 1x4 

E controlled variable matrix, 1x2 

e(t) error vector, 4x1 

F feedback gain matrix, 2x4 

G noise gain matrix, 4x1 

g standard free-fall acceleration, 9.80665 m/sec^ 

H matrix solution of equation (18) 

H(cj) frequency response function 

I identity matrix, 2x2 

i = 

J performance index 

K Kalman-Bucy filter gain matrix, 4x1 

ky radius of gyration about Y-axis, m 

L scale of turbulence, m 

J?y vane distance forward of airplane center of gravity, m 

M mass, kg 

Mq = ^oo^w^ m-N-sec^/rad^-kg-m^ 

Mq! = m-N/rad-kg-m2 

Mky2 
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M5e 


^6f 

Ni(t) 

No(t) 

P 

Q,R>S 

q(t) 

Ri 




S-w 

t 

u(t) 

V 

VI 


Vi 

vq 

Wg(t) 


= — — *^°° , m-N/rad-kg-m^ 

Mky^ 

= — ^ m-N/rad-kg-m^ 

Mky^ 

input white -noise process, rad/sec^ 
measurement white -noise process, rad 
normal acceleration, a 2 ;(t)/g 
matrix solution of equation (14), 4 x 4 

matrices defined by equation (11), 4x4, 2x2, 2x2, respectively 

pitch rate, rad/sec 

dynamic pressure, N/m^ 

diagonal weighting matrix, 2x2 

diagonal elements of matrix Rj 

wing area, m^ 

time, sec 

input vector, 2x1 

airspeed, m/sec 

matrix defined following equation (18), 4x4 
constant intensity of input noise process, Nj(t), m^/sec^ 
constant intensity of measurement noise process, NQ(t), rad^ 
vertical component of gust velocity (positive upward), m/sec 
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X,Y,Z 

WT 

Zq, 

Z6f 

«(t) 

6e(t) 

6f(t) 

6v(t) 

v{t) 


^ 1 , 2 , 3, 4 
?(t) 


body axes 

state vector, 4x1 

estimate of x(t), 4x1 

^oo^W 


da MV 

f£z 

36q MV 

^oo^W 


, N -sec/rad -kg-m 
, N-sec/rad-kg-m 
, N-sec/rad-kg-m 


85f MV 
angle of attack, rad 

scalar multiplier 

elevator deflection angle, rad 

flap deflection angle, rad 

vane deflection angle, rad 

variable appearing in equation (4), rad/sec 

variable of polynomial, rad/sec 

zeros of polynomial, rad/sec 

variable appearing in equation (4), rad 

variance of Wg, m^/sec^ 

rms normal acceleration 

rms pitch rate, rad/sec 


'a 


rms angle of attack, rad 
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^5e 

^6f 

^5v 

T 

•^>Wg(t)(w) 

CD 

I I 

[] 

E() 

Subscripts; 

e 

f 

g 

I 

max 

O 


rms elevator deflection angle, rad 
rms flap deflection angle, rad 
rms vane deflection angle, rad 
time increment, sec 

power spectral density function for input noise process (=vi), m2/sec5 

power spectral density function for measurement noise 
process (=Vq), rad^ 

power spectral density function for Wg, rad^/sec 

circular frequency, rad/sec 

absolute value of a quantity 

rectangular matrix 

expected value 

elevator 

flap 

gust 

input quantity 
maximum value 
measurement quantity 
vane 
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w 


wing 


free -stream value 
Superscripts: 

T transpose 

— vector 
" estimate 


transformed quantity 

Dot over a quantity denotes derivative with respect to time. 

THEORY 


Optimum design considerations of a gust alleviation system for an airplane cruising 
in turbulent air are studied in this paper. The elevator and a wing flap (that is, a trailing- 
edge control capable of positive and negative deflection angles) were used as controls to 
produce the force and moment necessary for reduction of normal acceleration caused by 
gusts. A measured vane deflection angle actuated the elevator and flap. A quadratic per- 
formance index should be minimized. Since atmospheric turbulence is random, the per- 
formance index was the variance of the normal acceleration and the control deflection 
angles. In addition, the airplane with the alleviation system in operation must be stable. 
Stochastic control theory was applied to obtain the desired alleviation system. 


Mathematical Model of Airplane Motion 

Linear equations of longitudinal motion that approximate the short period mode were 
used in this study. The frame of reference for the airplane motion is the system of body 
axes illustrated in figure 1 . The equations are similar to those given in reference 3 . The 
resulting differential equations follow: 


i(t) = Z^ait) + q(t) + i Zo;Wg(t) + Z5^6e(t) + Zgj6f(t) 

q(t) = MQ,Q!(t) + Mqq(t) + ^ MQ,Wg(t) + M6g8e(t) + M6j6f(t) 


( 1 ) 


J 
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For simplicity, the controls 6e(t) and 6f(t) were assumed to respond instantaneously 
when actuated. The influence of turbulence on the airplane motion was the vertical gust 
velocity Wg(t). Since Wg(t) is a random variable, it is only known statistically. Its 
mathematical representation as a function of time is discussed in a later subsection. 

Vane Deflection Angle 

The vane, located ahead of the airplane wing, is shown in figure 2. It is positioned 
so that it is not appreciably influenced by the flow field generated by the fuselage and wing. 
The sign convention used is also shown in figure 2. The angular displacement of the vane 
is represented by the following equation: 

6v(t) = -a(t) + ^ q(t) - 1 Wg(t) + No(t) (2) 

The response of the vane was considered to be instantaneous. 

In equation (2), the term NQ(t) has been included to account for the noise (uncer- 
tainty) in the measurement of 6y(t). The noise NQ(t) is assumed to be a stationary 
white -noise process having intensity vq and zero mean. Additional properties of 
No(t) necessary for this study are the covariance function 

E|No(t) No(rj> = Vo5(t - r) 

where 6(t - r) is a Dirac -delta function; and the power spectral density function 
for No(t) 


^No(t)(") = VQ 

which is constant for all frequencies w. Further discussion of white -noise processes 
may be found in references 6 and 7. 

Normal Acceleration 

Normal acceleration at the airplane center of gravity is 

az(t) = -v[q(t) - (i(t)] 

Using the first of equations (1), the normal acceleration is expressed by the equation, 
azW = vjzQ,a:(t) + ^ ZQ,Wg(t) + Z5g5e(t) + Z5f3f(t) 
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For convenience, the normal acceleration a 2 (t) is nondimensionalized by dividing by the 
acceleration of gravity g and is called normal acceleration n 2 (t) where 


nz(t) = 


az(t) 

g 


and 


Zo!a{t) + I ZQ,Wg(t) + X Z5^6e(t) + ^ 


(3) 


Mathematical Model of Atmospheric Turbulence 


Atmospheric turbulence is assumed to have the following characteristics: 
(a) It is one -dimensional (vertical gusts); 


(b) The vertical gust velocity Wg(t), a random variable, has a normal distribution 

function and a zero mean; 

(c) Its power spectral density may be adequately represented by the Dryden function: 


^Wg(^) = 


Low, 


2 1 + 


3L^ 

y2 


CO^ 


V2 


(d) It is stationary, homogeneous, and isotropic. 

The random variable Wg(t) can be mathematically represented (suggested in ref. 5) 
as the output of a linear system which has a stationary white -noise process Nj(t) as its 
input. This representation has been derived in appendix A and is given by the following 
equations: 


|(t) = T](t) 

m = -^ ?(t) - ^ vit) + Ni(t) 

Li 


Wg(t) = 4(t) + ?7(t) 


where the intensity of Nj(t) is vj 


o-Wg2v3 

l3 


(4) 


(5) 


9 


statement of Problem 


The formulation of the problem is obtained by manipulating equations (1) to equa- 
tion (5), Equation (5) is substituted into equations (1), (2), and (3). The resulting equa- 
tions are combined with equation (4), and matrix notation is used to obtain the following 
equations: 

• ^ ^ 

x(t) = x(t) + [bJ u(t) + [g] Ni(t) (6) 

6v(t) =[cJ^ + No(t) (7) 

nz(t) =[d]^+[e]^ (8) 


where 


x(t) = 


a(t) 
q(t) 

m 

is the system state vector; 


u(t) = 


6 e(t) 

6f(t) 


is the input vector; 




"Q! 


M 


a 


1.0 


M. 


is the system matrix; 



L^/3 „ 



a 


0 1.0 

l2 ' L 
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Z6e 



Mse 

0 


0 


is the output matrix; 

ro‘ 


0 



ZSf 

MSf 

0 

0 


is the system noise gain matrix; 


i c] = 

-1.0 

Iv 

-1.0 

L\/3 

J 


V 

V 

y2 

is the measurement matrix; 

and 



[d] = 

ISl 

1 ^1“ 1 

0 

- Zfl, 

g “ 

\/3 L 
gV 

and 





[e] = 

V 7 
g Z6e 




are the controlled variable matrices. 




Z 


a 


Equatio n (6) is a system of linear, constant coefficient, differential equations with 
the controls u(t) and white noise Nj(t) as inputs. The vane deflection angle 6y(t) 
given by equation (7) is assumed to be a measured quantity that is corrupted by white 
noise No(t). Normal acceleration n^Ct) given by equation (8) is the variable to be 
controlled. Since the input Nj(t) is assumed to be a Gaussian random variable, the 
system variables x(t), 6y(t), and n^Ct) are also Gaussian random variables. The 
problem is to determine a control functional u(t) - f(6y(t)j which minimizes the vari- 
ances of normal acceleration and the required control. That is, the quadratic perfor- 
mance index to be minimized is 


J = E 


nz^(t) 


+ u(t) 



where 



is a 2 X 2 weighting matrix which is to be specified. 


(9) 
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Thus, 


( 10 ) 


( 11 ) 


Equation (10) may be put into a form more suitable for calculations by applying the trans- 
formations given in appendix B. The resulting equivalent form for equation (10) is the 
following: 

J = E|k^'r[9]j^ + u(t)T|R]u(t)'| 

The solution to the problem depends on a knowledge of both the intensity of the mea- 
surement noise Vq and the value of the elements of the weighting matrix |R;^ . A per- 
fect measuring instrument is one for which vq = 0, and an undesirable instrument is one 
for which vq is greater than the quantity to be measured. The matrix |r^ must be 
positive definite and symmetric; that is, |rJ > 0 and = f^l]> respectively. 

Solution of the Problem 

The solution of the preceding problem has two interacting parts. In one part, the 
optimum regulator, a feedback gain matrix is calculated that gives u(t) as a linear 
function of an estimate of the state vector. This estimate is denoted x(t). The second 
part of the solution, the Kalman-Bucy filter, gives the structure and gain of an estimator. 
The estimator gives x(t) from the measurement 6y(t) and the control u(t). 

First the control input u(t) is chosen according to the linear control law that mini- 
mizes the performance index. (See eq. (10).) This optimum control law is written 

^JW=-[F]I(t) (12) 


nz^(t) = x(t)'^[p]'^ + u^'^[eJ'^|^[d]^ + [Eju(t)J 
so that equation (9) may be rewritten as 

J = E|x(tVF[Q]x(t) + 2x(t)'F'’s]u(t) + u(t)^[R]u(t)| 

In equation (10) the matrices [q], fsjj and (rJ are defined as 

[q]=[d]T[d] 

[sJ = [dJT[e] ^ 

[R]= [r,] +[e]T[e] j 
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( 13 ) 


where the feedback gain matrix [f] is 
[f] = [r]-i([b]T[p]+[s]t) 

(See eqs. (11) and appendix B.) The matrix [p] is the positive, semidefinite solution of 
the following matrix Riccati equation. (See eqs. (11) and appendix B for matrix definitions.) 

imi * [SFH * K - [PlB][R]-'[Bp[p] , 0 (14) 


The equations needed to obtain the feedback gain matrix are independent of the stochastic 
characteristics of the problem. In addition, the measurement matrix [c] is not used in 
the equations. 

Second, the matrix |jK^ is chosen in the system 


x(t) = |]^x(t) + jjB]u(t) + 6y(t) - [c]x(t) 


(15) 


so that the variance of the error 


e(t) = x(t) - x(t) 


is minimized. 


Thus, 


E 




(16) 


is a minimum. For this part of the solution, it is assumed that the white-noise pro- 
cesses Nj(t) and NQ(t) are uncorrelated. 

The system described by equation (15) is a Kalman-Bucy filter for the system given 
by equations (6) and (7). The filter gain matrix 'jcJ which minimizes the variance of the 
error (see eq. (16)) is 


[Kj = f[Hj[cF (17) 

where is the positive, semidefinite solution of the following steady -state matrix 

Riccati equation 

[h1[a] + [aJT [h] + [Vi] - i [h][c7 [c][h] = 0 (18) 

where 


[Vi] = vi[G]rGjT 

Although the input matrix [b] is included in the Kalman-Bucy filter equations, it is not 
needed to determine the filter gain matrix [k]. 
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The feedback control input given by equation (12), the Kalman-Bucy filter given by 
equation (15), and the filter gain matrix given by equation (17) constitute the optimum gust 
alleviator. Figure 3 shows a diagram of the closed-loop system. 

The weighting matrix [r^ and the measurement noise intensity vq are arbi- 
trary parameters. These parameters are varied to evaluate the performance of the gust 
alleviator. 


SAMPLE APPLICATION OF THEORY 

The theory outlined in this paper was illustrated by applying it to the optimum 
design of a gust alleviation system for a STOL airplane. The airplane mass, dimensions, 
flight condition, and aerodynamic characteristics used in the calculations are presented 
in table I. The flight condition was for cruise at an airspeed of 109 m/sec and at an alti- 
tude of 3048 m. The scale of turbulence was chosen to be L = 305 m and the mean 
square gust intensity was chosen to be = 1 m^/sec^; the corresponding noise inten- 

sity was vj = 0.04559 m^/secS. 

A wide range of values for the weighting matrix jjlj and the measurement noise 
intensity vq were used to calculate the following quantities: rms angle of attack ctq,, 

rms pitch rate Oq, rms elevator deflection angle ogg, rms flap deflection angle 05 ^, 
rms normal acceleration <^nz> estimated rms angle of attack ctq,, estimated pitch 
rate Oq, estimated rms gust velocity gain matrices [f] and [k], poles of the 

closed -loop system, and percent alleviation r. The rms of the variables a, q, n^, 

Wg, 5q, and 6 ^ were directly proportional to the rms gust velocity. Therefore, 
the results are valid for values of the rms gust velocity other than unity. A computer 
subroutine package entitled ORACLS (see ref. 7) was used for the necessary computations. 

Weighting Matrix Selection 

Selection of the weighting matrix |R]^ reflected the desire of the designer to keep 
the rms normal acceleration as near zero as possible with rms control deflections as 
small as possible. In most practical applications , jR^j was chosen to be diagonal. This 
choice permitted individually weighting 6 g and relative to normal acceleration. As 
an initial selection for this example, the elements rj and T2 were defined as the square 
of the ratio of the normal acceleration to the maximum deflection angle of the specific con- 
trol, Thus, 



14 



and 


■■ 2 = 

,max/ 

In these expressions, the acceleration n^ results whenever the maximum deflection 
angle of either 6^ or is the only input to the airplane. Adjustments can be made 
in rj and V 2 to satisfy the designer’s specifications as nearly as possible. The 
values of r j and V 2 computed for the airplane of this study are: 

r^ = (22.07)2/rad2 

V 2 = (16.60)^/rad^ 

It is not unreasonable to choose rj = r2 so that RjJ may be written as 
[Ri] = ^[l] 

where j3 is a scalar and TlJ is the identity matrix. Initially, the value of /3 was 
jS = (20)2/rad2. However, calculations were made for a wide range of values of /3 and 
for two values of measurement noise intensity vq. Results are presented in figure 4. 
The plots of and asj in figure 4 are of interest. The plot of Uq is pre- 

sented for convenience. It is evident that the value /3 = 400/rad^ was large enough to 
limit the control movements to small rms deflection angles and, consequently, rms nor- 
mal acceleration was relatively large. A much smaller value of j3 (jS < 10/rad2) would 
be more desirable for this problem since the rms control deflection angles would not be 
excessive for large gusts, and normal acceleration could be reduced substantially. For 
very large values of /3, the penalty for using the controls is so great that the alleviation 
system cannot produce large reductions in the rms normal acceleration. A value of 
|3 = S/rad^ was arbitrarily chosen for the system. 


Design Performance 

The performance of the gust alleviation system was measured by the reduction 
in rms normal acceleration. The unalleviated airplane, at the flight condition of this 
study, had an rms normal acceleration of 0.07928. Percent alleviation r may be cal- 
culated by using the following formula: 


0.07928 - CTnz 
^ " 0.07928 


100 


(19) 


where is the rms normal acceleration. 
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The Kalman-Bucy filter provided the estimate of the state x(t) needed in the 
optimum control law given by equation (12). Calculation of the Kalman-Bucy filter gain 
matrix |j^ depended on the value of the intensity vq of the measurement noise. (See 
eq. (17).) The gain matrix | i^ has been calculated for a wide range of values of the 
measurement noise intensity Vq with ^ = 3/rad2. The corresponding percent allevi- 
ation r is presented in figure 5 as a function of the measurement noise intensity. The 
ratio of the measurement noise intensity to the rms vane deflection angle vq^Oq^ is 
also shown in figure 5. 

Reference 1 showed that perfect alleviation theoretically can be achieved, with finite 
gains, for an airplane using a vane sensor and two controls and 6f). Perfect allevia- 
tion was not possible, under the assumptions of the present theory, when noise was super- 
imposed on the vane angle measurement. The measurement noise was propagated through 
the feedback loop and influenced the percent alleviation. However, the Kalman-Bucy filter 
reduced the effect of the measurement noise since the variance of the error vector e(t) 
was minimized. (See eq. (16).) Figure 5 shows that alleviation was 92 percent for a very 
small value of the measurement noise intensity (vq = 4.56 x 10 ”7 rad and 
vq/o^Sv ~ ^ lO'S). The corresponding filter gain is 


-5.4964 

11.5172 

[K] = 

7.1569 


( 20 ) 


-316.0585 


and the required rms control deflection angles are 
= 0.003485 rad^ 

> 

CTgj = 0.008179 rad i 


( 21 ) 


A more realistic instrument for measuring the vane deflection angle could have a 
noise intensity vq = 4.56 x 10“"^ rad or a ratio = 0.036. This ratio means 

that the noise is about 3.6 percent of the vane deflection angle. The optimum gust allevi- 
ation system then provides about 63 percent alleviation of normal acceleration. Table II 
gives a complete set of calculated characteristics for this optimum gust alleviation sys- 
tem* The Kalman-Bucy filter produced very good estimates of the state variables which 
are reflected by the good percent alleviation. The gains and are moderate 
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and the required rms control deflection angles are smaller than required for better 
alleviation. (See eq. (21).) 

If the possibility of the existence of ±20 percent tolerance in the rms measurement 
noise intensity is considered, this characteristic of the instrument could result in reduced 
performance of the alleviation system having the gains []f] and jjc] fixed. However, 
the calculated reduction in percent alleviation for this "off -design” condition was 
negligible. 

Table II includes the zeros of the polynomials in X given by the following 
equations: 


det |xT] - = 0 (22) 

and 

det |*X[lJ - >] + Jc;[c]j = 0 (23) 

The zeros of equations (22) and (23) are the poles of the closed-loop systems. (See fig. 3.) 
The zeros of equation (22) are the optimum regulator poles; the zeros of equation (23) are 
the Kalman-Bucy filter poles. Since the real parts of all the zeros are negative, the 
closed -loop system is stable. This result is not surprising since the regulator can always 
be stabilized with an additional feedback gain matrix, and the Kalman-Bucy filter is auto- 
matically stable. 

It should be emphasized that only one flight condition has been considered. In an 
actual design study, consideration of several flight conditions is essential. Also, in the 
problem formulation, control servodynamics must be included for a more accurate model 
of the system. The effect on alleviation associated with servo lag times can be partially 
offset by proper location of the vane. 

CONCLUDING REMARBCS 

A study has been made of optimum design considerations of a gust alleviation sys- 
tem for aircraft flying in turbulent air. A vane sensor (with noise) was used to measure 
vertical gusts, and elevators and flaps were used to reduce the root-mean-square value 
of the normal accelerations associated with the aircraft response to gusts. Since tur- 
bulence has stochastic properties, stochastic control theory was used in the analysis. 

A performance -index function involving normal acceleration and control deflections with 
a weighting matrix was used. 

A short take-off and landing (STOL) airplane in flight through turbulent air was used 
as an example to illustrate the application of the analysis. Effects of varying the noise 
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characteristics of the vane sensor and the weighting matrix in the performance -index 
function were determined. 

Calculations were performed as required by stochastic control theory to obtain the 
root -mean-square response of the airplane to turbulence. The calculations showed that 
a weighting number appearing in the performance index should be less than 10. A value 
of 3 was selected for subsequent calculations. Normal acceleration could be reduced by 
92 percent when the intensity of the measurement noise had a very small value. However, 
the filter gains were large and large control deflection angles were required. Good allevi- 
ation was calculated when the intensity of the measurement noise was about 3.6 percent of 
the vane deflection angles. Normal acceleration was reduced by 63.2 percent for moder- 
ate values of the gains. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
January 22, 1976 
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APPENDIX A 


MODEL OF ATMOSPHERIC TURBULENCE 


Development of a mathematical representation for the one -dimensional vertical 
gust having a power spectral density represented by the Dryden function was desired. 

The Dryden function has a scale length L and a mean square gust intensity Owg^. The 
desired representation can be developed as follows: 

In a linear system given by the following equations 


iit) = 7](t) I 

??(t) = -aQ4(t) - aj7?(t) + Nj(t)J 
the function Nj(t) is the input. An output is defined as 


(Al) 


Wg(t) = rj4(t) + 72 v(t) (A 2 ) 

It is desired to determine yj, 72 , aQ, and a^ to obtain the desired form of the 
gust power spectral density function. If Nj(t) is a stationary white -noise process, then 
?(t), ’iW, and Wg(t) are stochastic variables. Further, Nj(t) is assumed to be a sta- 
tionary white -noise process having the following characteristics: 

The mean of Nj(t) is 
E|Ni(t)j = 0 

and the convariance function is 

E |Ni(t) Ni(r)J = vi6(t - r) 

where Vj is the intensity of the process and 6 (t - t) is a Dirac -delta function. 

The power spectral density function is 


^Ni(t)(w) = Vi 

which is constant for all frequencies w. 

The power spectral density function of the output can be calculated by the following 
relationship: 
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^Wg(t)(w) = iH(o)) ^%j(t)(w) 


or 


# 


W| 




H(o-) 


n 


(A3) 


where |h(co)| is the absolute value of the complex frequency response function H(w) of 
the output Wg(t), The frequency response function of Wg(t) is (from eqs. (Al) and (A2)) 


and 


H(o.) = 


H(w) 


- 0)2 


+ iwaj^ 

yj2 + 


2o;2 


w'i + - 2ao)w2 + 

Consequently, by using equation (A3) 


$ 


w. 


r(t)(“) = 


[ri‘ 


y 2 


2^,2' 


co^lv 


I 


+ (a-1^ - 2ao)o)2 + ag2 


(A4) 


Since Wg(t) may be considered the vertical component of gust velocity in atmospheric 
turbulence , evaluation of the parameters y^, y 2 , ag, and a^ may be done in the fol- 
lowing way. Atmospheric turbulence may be assumed to be a one -dimensional, isotropic, 
stationary, and homogeneous random process. Also, the Dryden power spectral density 
function 


4> 


w, 


,(t)(^) = 


^Wg L 




1 + 


V2 


(A5) 


may be considered to represent atmospheric turbulence. The functions given by equa- 
tions (A4) and (A5) are equated and the constants and coefficients of co are equated; 
the following values for the unknown parameters result. If the constant is set 
equal to 1.0, then 


yi 


= 1 


L{3 

y 2 = — 


^0 - 


l2 
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Vi = 


Owg^vS 

l3 


Note that > 0 and a^^ > 0 are used to assure stability of the system of equa- 
tions (Al). These parameters are substituted into equations (Al) and (A2) and the follow- 
ing model of the assumed atmospheric turbulence is obtained: 


and 


?(t) = vit) 

2 / 
’7(t) = ^(t) - ^ Vit) + Ni(t) 


Wg(t) = |(t) + ^ 7]{t) 


(A6) 


(A7) 


One -dimensional atmospheric turbulence has been modeled as the output of a linear 
system with a white -noise process as the input. The power spectral density function for 
the vertical gust velocity is the Dryden function. 
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APPENDIX B 


SYSTEM EQUIVALENCE 


In optimum control problems, the quadratic performance index often contains cross 
products of the states and the controls. A transformation is presented here to show that 
such a performance index is equivalent to another one which does not contain the cross - 
product terms. The dynamics of a plant are assumed to be represented by the following 
constant linear differential equation: 

^(5 = [a]1^+[b]^ (B1) 

where x(t) is the state vector, u(t) is the input vector, [Xj is the system matrix, 
and [^b]] is the input matrix. A well-known result (see ref. 5 ) of optimum control 
theory is that the feedback control law given by 

of) = - [r] ■ ^ [b]T[p];^ (B2) 


minimizes a performance index given by 

J = E^^TgJ^ ^ (33) 

and the symbol E{ } indicates the expected value when x and u are random 
variables. 

The symmetric matrix [pj in equation (B2) is the positive -definite solution of 
the following matrix Riccati equation: 

0 , -[p][X] - [S]T[p] - [g] +[p][b][rJ-"[b]’'[p] 

Equations (Bl) to (B4) are called System I. 

The following transformations are applied to System I: 

^ ^ - [R]"l[s]'^x(t)' 

[A] = [A].[Bj[R]-l[sf (B5) 

[Q] = [Q]t[s][R]-l[s]T 
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and the result is System n, given by the following formulas: 


+ [b]uW 

= -I>T^[[bFW +[sF]^ 

j = + 2 j^[s]^[w +‘JIw'^[h]^ 

The corresponding matrix Riccati equation results: 


0 = -[P][A] - [a]T[p] - [Q] . ![p][b] . [stRj-lfCBf [P] . [S]T 


(B6) 

(B7) 

(B8) 

(B9) 
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TABLE I.- AIRPLANE MASS, DIMENSIONS, FLIGHT CONDITION, 

AND AERODYNAMIC CHARACTERISTICS 

Mass, M, kg 5670 

Wing area, S-^^,, 39 

Mean aerodynamic chord, c, m 1.981 

Radius of gyration about Y-axis, ky, m 2.572 

Vane distance ahead of center of gravity, £y, m 2.972 

True airspeed, V, m/sec 109 

Altitude, m 3048 

Dynamic pressure, q^o, N/m^ 5364 

Zq,, N-sec/rad-kg-m -1.969 

^6e> N-sec/rad-kg-m -0.156 

Zgj, N-sec/rad-kg-m -0.746 

Mq,, m-N/rad-kg-m^ -14.597 

Mq, m-N-sec^/rad^-kg-m^ -2.095 

Mgg, m-N/rad-kg-m2 -20.042 

Msf, m-N/rad-kg-m^ 8.672 
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TABLE n.- CHARACTERISTICS OF AN OPTIMUM GUST ALLEVIATION SYSTEM 


/3 per rad2 . 



3 

Vq» rad 



4.56 X 10-'^ 

vo/<^6v 



0.036 

rad 




Oq,, rad 



0.007375 

Oq, rad/sec 



0.01527 

dq, rad/sec 



0.01522 

Owg, m/sec 



1.0 

owg, m/sec 



0.8880 

<^6e’ 



0.003268 

CT6f, rad 



0.007643 

%z 



0.02914 

r, percent 



63.2 

[-1.0405 

-0.2920 

-0.8337 

-4.2172~ 

li 

1 1 




|_2.7328 

0.0611 

2.6892 

13.0734_ 


-4.6441 


12.2582 

[K] = 

4.3912 

-9.4190 

Zeros of det [xi - A + BF^ = 0: 
Xi = -0.3573 
X2 = -0.3572 
X3^4 = -4.2838 ± 16.4486 

Zeros of det [xi - A + KC] = 0: 
Xj = -4.8199 
X2 = -2.5355 
X3 = -1.0100 
X4 = -1.9362 
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Figure 1.- Axis system. 


to 
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00 



Figure 2 


Vane location and sign convention. 












Basic airplane 



4. 56 X 10'^ rad 
a 28 X I0'2 rad 


I 



2 

Weighting number, (3, 1/rad 


(a) Root -mean-square (rms) normal acceleration and pitch rate 
Figure 4.- Airplane and gust alleviator response to turbulence. 



4. 56 X 10 rad 
^ 28 X I0’2 rad 



(b) Root -mean-square (rms) control deflection angles. 
Figure 4.- Concluded. 
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